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ABSTRACT 

We present the results of a study of the two-point correlation function for 
a sample of field galaxies taken from the CNOC cluster survey. The sample 
consists of 144 galaxies within a contiguous region of space subtending 225 
square arcminutes. The objects have r-band magnitudes 17.0 < r < 21.7 and 
redshifts 0.21 < z < 0.53. The median redshift of the sample is z = 0.36. 
The real space correlation function is found to be consistent with a power law 
£( r ) = ( r / ro )- L7 with r = 2.1±g;| fc^Mpc (fi„ = 1), or r = 2.5±g:| fc^Mpc 
(Qq = 0.2). Uncertainties are estimated using the bias-corrected bootstrap 
resampling method, with 300 resamplings. This low correlation length implies 
strong evolution since £~0.36 has occurred in either the correlation function 
or the luminosity function; if the observed correlation function is modeled as 
£{r,z) =£(r,0)(l + ^)-( 3+£ ) with £(r,0) = ( y 51 h . lMpc )~ 1 - 7 , then e = 0.8^;°. 
Comparison of the redshift space and real space correlation function indicates 
that the one-dimensional pairwise peculiar velocity dispersion a at z ~ 0.36 is 
weakly inconsistent with 770 kms -1 , the value predicted by the Cosmic Virial 
Theorem if Q = 1. The observed correlation function is, however, consistent 
with a = 400 kms" 1 , the value expected if Qq — 0.2. 
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1. Introduction 

The study of the statistics of galaxy clustering has yielded important information 
about the large-scale structure of the universe and about the environment of galaxies. One 
of the most useful statistics employed in this study is the two-point correlation function. 
This statistic quantifies the clustering of galaxies, and is directly related to the power 
spectrum of density fluctuations in the galaxy distribution. Determining the evolution of 
the correlation function is therefore essential for an understanding of cosmological structure 
formation. 

The two-point correlation function has been extensively studied at low redshifts; 
[Efstathiou (1995")1 gives a summary of redshift surveys that have been used for correlation 



analysis. Observations of the correlation function at the present epoch indicate that it is 
well described by a power law 



at scales r ^ 10 ft, -1 Mpc. Here, and throughout this paper, all separations are given in 
physical, as opposed to comoving, coordinates unless otherwise stated; the Hubble parameter 
is taken to be Hq = lOO/i kms~ 1 Mpc~ 1 . Typical results from optical surveys of nearby 
galaxies include r$ = 5.4 ± 0.3 /i _1 Mpc, 7 = 1.77 ± 0.04 from the Center for Astrophysics 
(CfA) survey (Davis & Peebles 1983) and r = 5.1 ± 0.2 h^Mpc, 7 = 1.71 ± 0.05 from 
the Stromlo-APM survey (Loveday et al. 1995). Values for r$ and 7 for local samples are 
typically taken to be 5.4 /i _1 Mpc and 1.8, respectively, based on the CfA observations. 
For the purposes of comparing our results to those from local samples, however, we adopt 
r = 5.1 /i -1 Mpc and 7 = 1.7, principally because of the fainter magnitude limits used in 
the Stromlo-APM survey. 

The principal aim of this investigation is to examine the evolution of the correlation 
function. A useful empirical model for this evolution, introduced by [Koo fc Szalay (1984)| is 

£(r,z)=£(r,0)(l + ^ r (3+£) . (2) 

In this model, e = 7 — 3 corresponds to clustering fixed in comoving coordinates, as seen in 
biased Cold Dark Matter (CDM) simulations (Carlberg 1991). If e — 0, clustering is stable 
in physical coordinates. |Colin fc Carlberg (1996 )| find s ~ for an Open CDM (OCDM) 
initial power spectrum with Qq = 0.2. 

To determine s, one must compute the correlation function from data from an earlier 
epoch. The most common method is to compute the angular correlation function from the 
angular galaxy distribution (e.g. Koo & Szalay 1984, Efstathiou et al. 1991 and Infante 
& Pritchet 1995). This method does not require that the distances to the objects in the 
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sample be known; rather, the clustering is observed as a two-dimensional projection of the 
three-dimensional clustering over a wide range of redshifts. In order to estimate the amount 
of clustering evolution present in the data, therefore, models for both the present-day 
correlation function and the redshift distribution of objects in the sample must be employed. 
Both Efstathiou et al. (1991) and Infante & Pritchet (1995) find e > if fi = 1 and there 
is no evolution in the luminosity function. 

In order to produce a reliable estimate of the spatial two-point correlation function, 
redshifts for a large number of objects must be obtained. Pencil-beam surveys using 
multi-object spectroscopy (MOS), such as the Autofib survey (Cole et al. 1994) and the 
Canada-France Redshift Survey (CFRS) (Le Fevre et al. 1996) are the only surveys to 
date which contain a sufficient number of objects at intermediate redshifts. Le Fevre et al. 
(1996) find £~0 — 2 in the CFRS survey, consistent with the evolution observed in angular 
surveys. 

In principle, data from pencil-beam surveys may be used to estimate the redshift 
space correlation function. The required velocity accuracy is, however, quite high, pavis] 



fc Peebles (1983)] find the one-dimensional pairwise peculiar velocity dispersion in the CfA 



survey to be a = 340 ± 40 kms -1 at a separation of 1 h~ 1 Mpc. If the velocity errors in a 
survey are comparable to u j computation of a reliable redshift space correlation function 
from that survey is impossible. Cole et al. (1994) find a comoving correlation length 
tq = 6.5 ± 0.4 h~ l Mpc based on the redshift space correlation function for the Autofib 
survey. The median redshift of this survey is z = 0.16; taking r (0.16) = 6 ' 5 /i.i6 h~ 1 Mpc 
and r (0) = 5.1 h~ 1 Mpc in equation (H) yields e ~ —4, at odds with previous angular 
correlation results. 

We present here estimates of the real space and redshift space two-point correlation 
function for a sample of field galaxies taken from the Canadian Network for Observational 
Cosmology (CNOC) cluster survey. The data are described in Section |2|. The techniques 
used for computing the real space two-point correlation function and its evolution, along 
with the results for this survey, are presented in Section 0. Also described there is 
our investigation of the redshift space two-point correlation function and the heretofore 
unknown pairwise peculiar velocity dispersion at z ~ 0.36. Our conclusions are given in 
Section f|. 



-4- 



2. Data 

The data used here are taken from the CNOC cluster survey. The observational 
procedures and data reduction are described in detail in |Yee, Ellingson fc Carlberg (1996) 



only the relevant features of the survey are discussed here. The data were obtained using 
the multi-object spectrograph (MOS) at the Canada- France-Hawaii Telescope (CFHT). A 
band-limiting filter was used with the spectrograph to reduce the length of each spectrum, 
further increasing the multiplexing rate. Since the CNOC survey is a cluster redshift 
survey, each field was chosen to contain a cluster at its center - most of the data is 
therefore unsuitable for investigating the clustering of field galaxies. One field in the survey, 
MS1512+36, however, is well suited to a correlation analysis of field galaxies, since the 
cluster is quite poor. Only data from this field are analyzed here; the complete data set is 
described and presented in |Abraham et al. (1996)| . 

Pencil-beam surveys present several problems related to the fairness of the sample (a 
sample is said to be fair if the structure contained within it is representative of the global 
average). The first difficulty relates to the beam width; a typical beam diameter of 10' 
yields a width of~ 1.8 /i _1 Mpc at z = l /z (VLq = 1), the same order of magnitude as the 
expected correlation length r at that redshift, given any reasonable amount of evolution in 
the correlation function. 

Thus, the distribution of objects in a beam may be dominated by a single large 
density inhomogeneity, biasing the estimation of £. One possible solution to this problem 
is to calculate the correlation function from data taken from many pencil-beams scattered 
randomly throughout the sky. This approach does not make optimal use of the data, 
however; it is preferable to place the beams so that the beam-beam separation is comparable 
to the beam width, thus increasing the number of pairs of objects with separations ~ tq. 
The MS1512+36 data are from a mosaic of three fields from the survey, with a total angular 
size of 27' x 8' (222 square arcminutes). Despite the relatively large width of the sample 
(4.9 /i -1 Mpc at z — 0.36 if Qq — 1), there remains a large overdensity spanning the width 
of the field. The effect of this inhomogeneity is discussed below. 

A second, more subtle, problem with MOS data is that of selection effects. Magnitude 
selection may bias the sample towards bright objects, which will lead to erroneous results 
if bright objects cluster differently from faint objects. More importantly, MOS produces a 
lower limit on the separation of objects for which spectra may be observed. Once one object 
is designated to be observed through a slit on a given mask, the placement of the spectra on 
the detector precludes designing another slit closely above or below the first. This results 
in high density regions being sampled less completely than low density regions, thereby 
reducing the observed correlation. Each of the three fields composing the MS1512+36 data 
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was observed with 2 different MOS masks, with a higher priority given in the second mask 
to completing observations of closely spaced pairs (see Yee et al. 1996 for details). This 
somewhat reduces the amount of geometric selection. 

To correct for the magnitude and geometric selection, the magnitude weight w m and 
local magnitude weight W[ m are calculated from the data for each object in the sample. For 
a given object with apparent magnitude in the bin (m, m + Am) and observed redshift z, 
w m is proportional to the fraction of objects anywhere in the sample in the same magnitude 
bin which have observed redshifts. The local magnitude weight wi m for an object is 
proportional to the fraction of objects in a circle with radius 120" about the first object 
in the same magnitude bin which have observed redshifts. Also defined for each object 
is the geometric weight w xy = wi m /w m , which is related to the number of nearby objects 
at any magnitude which have observed redshifts. A detailed explanation of the weighting 
procedure is given in Yee et al. (1996); a test of the extent to which these weights correct 



for the sampling nonuniformities is described in Section 3.1 



Three subsamples of the MS1512+36 data are created; the photometric, redshift and 
field samples. The photometric sample, used for computing the angular correlation function, 
consists of the 404 objects with g — r colours, with Gunn r-band magnitude in the range 
17.0 < r < 21.7. The upper limit is chosen so that the magnitude weight of every object in 
the redshift sample is less than 5; the lower limit is employed since the masks were designed 
to exclude objects much brighter then the brightest cluster galaxy, which has r = 18.45. 

The redshift sample consists of all objects in the photometric sample which have 
identified redshifts in the range 0.21 < z < 0.53. These limits are chosen so that the spectral 
features used to identify emission line objects lie within the optimal response region of the 
filter. This sample is shown in Figure |I| The use of band-limiting filters for the CNOC 
survey results in strong redshift selection effects in the sample, which are easily understood 
by noting the visibility of strong spectral features in our limited spectral window. In this 
analysis we have extended the lower z limit from 0.27 in |Yee, Ellingson fc Carlberg (1996) 



to 0.21. The original higher limit was based on the detection of [O II] A3727 at the blue end 
of the spectrum. However, the wavelength limits for the filter used are such that as [O II] 
A3727 disappears at the blue end, the [O III] AA4959, 5007 lines come in on the red end. 
From our sample, we have found that whenever the [O III] AA4959, 5007 and [O II] A3727 
lines are both within the spectral range, they are always detected simultaneously. Hence we 
can safely use the [O III] AA4959, 5007 lines to extend our lower redshift limit; no significant 
selection bias as a function of spectral type is seen in the redshift sample (Yee et al. 1996). 

The field sample is constructed from the redshift sample by removing the 38 objects 
that have redshifts 0.3656 < z < 0.3796. This redshift range corresponds to a velocity range 
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of ±2100 kms" 1 , or six times the velocity dispersion of the cluster, at z — 0.3727, the center 
of the cluster (Carlberg et al. 1996). The field sample therefore almost certainly excludes 
all cluster members, thus removing the bias towards high density regions present in the 
redshift sample. The field sample contains 144 objects, with a median redshift of z — 0.36. 

As can be seen in Figure 1, the cluster appears to be embedded in an overdensity 
extending across the field. In order to test the sensitivity of our results to the presence 
of this structure, we have computed the correlation functions for the field sample with all 
objects with 0.354 < z < 0.390 removed. The results are consistent with those computed 
from the entire field sample, although, as there are only 118 objects in the field sample 
outside the these extended redshift limits, the uncertainties are considerably larger. We 
conclude that the presence of this structure does not unduly influence our results. 



3. Analysis 
3.1. Estimating the Correlation Function 

The two-point correlation function £ is defined by (Peebles, 1980) 

5P( r ) = n(l + £(r))5V , (3) 

where 5P(r) is the probability of finding a second object in a volume SV with a physical 
separation r from a randomly chosen object, and n is the mean density of objects. For 
a finite sample of objects, and some (small) fixed separation difference Ar, £(r) may be 
estimated from equation (13) as 



DD(r) V 

1 + « r )^AFw' < 4 > 



Here, Njj is the number of objects in the sample, DD(r) is the number of ordered pairs of 
objects in the sample with separation between r and r+Ar, V is the volume of the sample, 
and AV(r) is the average volume of the set surrounding a first object in which a second 
may be found with separation between r and r + Ar from the first. 

The volumes AV in equation (|j) may be arbitrarily complicated; in practice, they 
are estimated using Monte-Carlo integration. A random data set containing Nr objects is 
generated within the volume of the original data set in a manner such that the random 
catalog is subject to the same selection criteria as the data. If the number of pairs of objects 
with separation between r and r+Ar, the first object belonging to the data set, the second 



- 7- 



to the random set, is DR(r), then (Davis & Peebles 1983) 

, , N R DD(r) 

The pair counts DD and -D.R in equation (|J) may be computed with arbitrary weights, 
so that DD(r) = J2i,j w i D ^ w j D ^ an d DR(r) = J2i.j w i D ^ w j R \ where the sums are taken 
over all data-data or data-random pairs of objects with separation between r and r + Ar, 
respectively, and w\ D ^ and w\ R ^ are the weights to be applied to z'th data and random 
object, respectively The object counts No and Nr in equation @ are replaced with the 
weighted object counts D = ] 
function is then estimated as 



weighted object counts D = Y^d=xw\ D ^ and R = J2i=i w i R \ respectively. The correlation 



, , RDD(r) 
l + ■ (6) 

To correct for the selection effects present in this sample, and for the incompleteness present 
in any magnitude-limited sample, we take the total weight for an object at redshift z with 
local magnitude weight w\ m to be w to t = Wi m /<p(z), where <f>(z) is the redshift selection 
function, defined as the fraction of objects at redshift z which lie within the apparent 
magnitude limits of the survey. These total weights are then used for computing the pair 
and object counts defined above. 

The selected random catalog is generated by first creating a uniform random catalog, 
by randomly distributing objects throughout the sample volume in such a way that the 
comoving density of objects is constant. Apparent magnitudes are then assigned to each 
object in the uniform random catalog using the process described below, and those objects 
lying outside the magnitude limits of the survey are discarded. Finally, the local magnitude 
selection is estimated at the position of each object random object, and objects failing the 
selection criteria are discarded. For small data sets, this method is preferable to smoothing 
the observed redshift distribution, since the redshift distribution is dominated by density 
inhomogeneities comparable to the width of the sample, making the random distribution 
obtained sensitive to the smoothing window used. 

Absolute magnitudes for the objects in the field sample are required to determine the 
luminosity function, which is used for generating the absolute magnitudes for the random 
objects. The K-corrections are obtained by interpolating from model K-corrections in the r 
and g bands as a function of redshift for non-evolving galaxies of 4 spectral types (E+SO, 
Sbc, Scd, and Im). The models are derived by convolving filter response functions with 
spectral energy distributions in Coleman, Wu, & Weedman (1980). These values are then 
corrected from the AB system to the standard Gunn system (Thuan & Gunn 1976). For 
each galaxy with redshift, a spectral classification is estimated by comparing the observed 
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g — r color with the model colors at the same redshift. The spectral classification, obtained 
via interpolation, is treated as a continuous variable between the 4 spectral types. From 
the spectral classification, the appropriate K-correction to the r magnitude is then derived 
using the models. 

The r-band absolute magnitudes for the random catalog are generated according to 
the luminosity function for the data. The luminosity function is modeled as a non-evolving 
Schechter function (Schechter, 1976) 

$(M r ) = 0.41b lO^lO^+^W-^exp (-10°' 4 ( M *-^) , (7) 

where $(M r ) is the comoving number density of objects with absolute magnitude M r per 
unit magnitude. The assumption of no evolution is reasonable, since the redshift range of 
the field sample is relatively small. The parameters M* and a are adjusted to fit the data 
using a least-squares fit. The luminosity function (after an appropriate renormalization, 
which removes the dependence on </>*) is then used as the probability distribution for the 
absolute magnitudes of the objects in the uniform random catalog. An intrinsic colour 
{g — r) is then chosen for each object, in a manner such that the distribution of intrinsic 
colours in the uniform random catalog is the same as that in the data catalog. The r-band 
absolute magnitude M r and intrinsic colour (g — r) for each object are then used to 
compute the K-corrections, using the method described in Section [2|. 

For each object in the uniform random catalog with apparent magnitude within the 
sample limits, magnitude and geometric weights w m and estimated by interpolating 

from the data weights. The magnitude weight interpolation is a straightforward one- 
dimensional interpolation. The two-dimensional geometric weight interpolation is performed 
by convolving the spatial map of geometric weights of each data object with a Gaussian 
with dispersion o = 1'. The values of the convolved weights at the position of the random 
object are then summed, yielding a weight for that random object. The local magnitude 
weight wim for the random object is than calculated from w m and w xy ; the object is 
discarded unless l/wi m is less than a randomly chosen number between and 1. The 
redshift selection function (p(z) is then computed from the luminosity function parameters 
M* and a and the sample magnitude limits; wi m and <p(z) are then combined to give a total 
weight wtot = Wi m /<p(z) for each remaining object in the selected random catalog. 

The redshift distribution for a random catalog generated in this manner is compared 
to the redshift distribution of the redshift sample in Figure |2|. The random distribution has 
been rescaled to have the same number of objects within the redshift sample redshift range 
as the data distribution in this figure. Note that the random distribution is much more 
uniform than the distribution that would be obtained by smoothing the data distribution, 
due to the presence of the large structure discussed earlier. 
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This method of generating the random catalog and the weighting procedure are tested 
by comparing the angular correlation functions of the photometric and redshift samples. 
The angular correlation function w is estimated in a manner similar to that given in 
equation (|S|); 

, m R DD{6) 
l+W ^^DDW)^ (8) 

where the pair counts are now taken over pairs with angular separation 9. If the model (|l|) 
for £ is correct, then (Peebles 1980) 



w{0) = A W -) , (9) 




where £ = 7 — 1. 

The random catalog used for computing w for the redshift sample is generated using 
the procedure described above, while the catalog for the photometric sample was generated 
uniformly; both catalogs contain 10 000 objects. Figure shows the results of a power-law 
fit to the angular correlation for the photometric sample, with 5 fixed to 0.7, and w for 
the redshift sample. The amplitudes in equation (|9|) are found to be A w = 1.3loa f° r the 
photometric sample, and A w = 1.2lto]o f° r the redshift sample. The consistency of these 
amplitudes indicates that the total weights adequately describe the selection effects present 
in the data. The photometric sample value of A w implies to(l') ~ 0.07, consistent with 
the value found by Infante & Pritchet (1995) for a sample of objects with F < 22. The 
uncertainties in all correlations and fitted parameters given here are the 68.3% confidence 
intervals computed using the bias-corrected bootstrap resampling method, described by 
Efron fc Tibshirani (1986"]]. Three hundred resamplings are used for each computation, the 



minimum recommended for this method. 



3.2. The Real Space Correlation Function 

Two different correlation functions are derivable from redshift data; the redshift space 
correlation function £(s) and the real space function £(r). If the separation of a pair of 
objects is computed directly from redshifts, it includes the line-of-sight component of the 
object's peculiar velocities relative to the Hubble flow. For the range of separations of 
interest here, the random internal motions of bound groups of objects dominate, elongating 
structures along the line of sight. This elongation reduces the observed correlation at small 
separations, so the power-law model (|l|) is expected to be valid only for the real space 
correlation function. 
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3.2.1. Method 

Although it is impossible to measure real space separations directly using only redshift 
data, it is possible to estimate the parameters of a model for the real space correlation 
function. This is accomplished by decomposing the redshift space separation of a pair of 
objects into components parallel and perpendicular to the line of sight to the pair. Since 
the redshift space distortions act only along the line of sight, functions of the perpendicular 
component must be independent of these perturbations. 

The decomposition is performed assuming the separations are small, so that the effects 
of curvature may be neglected. Thus, given two objects with redshifts z\ and z 2 and angular 
separation 12 , two vectors Xj (i = 1, 2) are formed, such that 



2c tt z - (2 - ftoXvT + TV - K 



1,2) (10) 



(A = 0), and 

Xi • X 2 = |xi| |x 2 | COS012 • (11) 

The comoving redshift space separation of the pair is then s ~ x 2 — xi, and the line of 
sight to the pair is x = | (xi + x 2 ). The components of the physical separation parallel and 
perpendicular to the line of sight are then 

* = 7t4^fi (12) 

(1 + z) |x| 

and 



|2 2 
S — Tt Z 



r " = \ + T~ ' (13) 

where z — | {z% + z 2 ). The definition of the correlation function (^) is then generalized 
to describe the probability in excess of random of finding an object with redshift space 
separation (r p , it) from a randomly chosen object. This is estimated as 

R DD(r v , it) , , 

where DD(r p ,7r) and DR(r p ,7t) are the weighted number of data-data and data-random 
ordered pairs with separations (r p ,7t), respectively, and D and R are the weighted object 



counts defined in section ETT 



Although £(r p ,7r) is affected by the redshift space distortions described earlier, the 
projected correlation function w p (r p ), defined by 



w P ( r P ) = / €(r P ,n)d7t (15) 
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is not. Thus 



w P ( r p) = € (\Ap + x2 ) dx ' ( 16 ) 

where the integral is over the real space correlation function. If the power law model ([!]) is 
employed, then (Davis & Peebles 1983) 



= ^Ym\ ro (%) ■ (17) 

The integral in equation (|l^) must be truncated at some n max for any real data set. 
Figure ^ shows the quantity 27rnJ p (l /i _1 Mpc, tt) vs tt for the field sample, where 

J P (r P , tt) = £ jf* £(r' p , 7r')r»' . (18) 

Given this definition, the quantity 27tnJ p (r p , tt) represents the mean number of objects in 
excess of random within the cylinder with radius r p and length 2tt centered on an object 
in the sample. This function is expected to increase with tt for small tt, approaching some 
limiting value. As can be seen, cutoffs less than ~ 3 /i _1 Mpc exclude real power in £(r p ,7r), 
while noise appears to be the primary contributor to the integral ( |TB"D for tt ^ 50 /i _1 Mpc. 
We adopt 15 /i _1 Mpc as the cutoff to be used in equation fllBl); no significant change in 



-s 



the derived correlation length is observed when the cutoff is varied between 5 h 1 Mpc and 
35 h^Mpc. 



3.2.2. Results 

The projected correlation function w p {r p ) is calculated for the field sample using 
equation (14), with the integral truncated at 15 /i -1 Mpc. The random catalog is generated 



using the method described in Section |3Tlt the catalog contains 200 000 objects with 
redshifts. The pair counts DD and DR and object counts D and R are calculated using the 
total weights defined in section |3.1| . The results are shown in Figure |5|; r is determined for 
Qo — 1 and 0.2 by fitting the model ([17]) to the data with 5 fixed to 0.7. The correlation 
length is found to be r = 2.llgj h^Mpc (Q Q = 1) or r = 2.5^ h^Mpc (Q Q = 0.2), 
considerably smaller than the values found locally. 

If it is assumed that the population in our sample at z ~ 0.36 will evolve to the 
population observed at z = in optical surveys, equation (H) may be used, with r (0) 
equal to the value observed in the local sample, to estimate e. Applying equation (fj), with 
r = 2.118:1 /i _1 Mpc at z = 0.36 and r = 5.1 h^Mpc at z = yields e = 1.8±?;§ (fi = 1); 
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using r = 2.5ioi h~ 1 Mpc gives e = 0.8lj;;3 (f2 = 0.2). Thus, the correlation is found to 
be increasing with redshift in physical coordinates, although the Qo = 0.2 estimate of e 
is consistent with no evolution. Correlation fixed in comoving coordinates (e = 3 — 7) is 
effectively ruled out. 

An alternate explanation for the low correlation length is that some of the objects 
in our sample are weakly clustered and become intrinsically faint at the present epoch 
(Efstathiou et al. 1991). That is, the evolution is in the luminosity function, not the 
correlation function. This is entirely feasible in our sample; if Tq — 5.1 /i _1 Mpc locally, 
and r = 2.5 /i -1 Mpc at z = 0.36, and e = 0, and the faint population is completely 
unclustered (£ = 0), then the fraction of objects in our sample which belong to this 
currently faint population is only ~ 0.15. If the faint population is clustered realistically, say 
ro = 3.8 /i -1 Mpc (as seen in the IRAS survey, Fisher et al. 1994a) for both the faint galaxy 
autocorrelation and faint-bright cross-correlation, then ~ 40% of the objects in our sample 
are required to be undetected in the present-epoch observations. Thus, we cannot rule out 
evolution in the luminosity function as the source of the reduction in r$ at intermediate 
redshift s. 



3.3. The Redshift Space Correlation Function 

As noted in section |T^, the redshift space correlation function is not expected to be a 
power-law due to random peculiar velocities. However, measurements of £(s) can be used, 
in conjunction with a model for the real space correlation function, to provide information 
about the velocity distribution of objects in the sample. This in turn yields information on 
the mean matter density. 



3.3.1. Method 

The real space correlation function £(r) is related to £(r p ,7r) by (Peebles 1980) 

1 + t(r p , tt) = J g(r, v) (1 + £(r)) a?v , (19) 

where g(r,v) is the distribution of relative peculiar pairwise velocities of pairs with 
separation r, and r 2 = r 2 , + (it — v z j H(z)) 2 , where v z is the component of v along the line 
of sight. To obtain a relation between the redshift space correlation function £(s) and the 
real space correlation function £(r), one integrates equation ( |j~9|) over a sphere of radius s; 

Air J" £(s')s' 2 ds' = 2tt J J* ^ " £(r)g(r,v)r r dr p dTrd 3 v , (20) 
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using the identity 2ir J* s ^ £(r p , -K)r p dr p dix = Atx J s £(s')s' 2 ds' . 

We employ here a simplified model for the pairwise peculiar velocity distribution in 
which g is independent of the separation r. Two dimensions of the velocity integral in 
equation (|20|) may therefore be performed immediately; we define the line-of-sight peculiar 



pairwise velocity distribution f{y z ) = Jf^ g(v x , v y , v z )dv x dv y . Equations (|20D then 
reduces to 

4tt / Z(s')s' 2 ds' = 2tt / / / i(r)f(v z )r p dr p dTidv z , (21) 

JO J-ooJ-sJO 

where the integral on the left is over the redshift space correlation function, while the 
integral on the right is over the real space correlation function. Differentiating equation 
(f2"T|) with respect to s gives the general relationship between the redshift space and real 
space correlation functions, under the assumption that g is independent of r; 

£(s) = l^JZfj {^ 2 -2{v z /H{z))it + {v z /H{z)) 2 ^ f(v z )dndv z . (22) 

The argument to the real space correlation function in equation ( p2]) is just the physical 
separation r, evaluated with r p = \/ s 2 — 7r 2 . If £(r) is modeled as a power law (0), then the 
integral over 7r in equation ( ^2] ) may be performed analytically, finally yielding 

£(s) = ^^(^o*- 1 (I* + v z /H(z)\ 2 ^ -\s- v z /H(z)\ 2 ^) f(y,)^ . (23) 

The simplified model for the velocity distribution used here takes / to be an exponential 
with zero mean, and dispersion independent of separation; 



= V 2^ 6XP ( 



(24) 



a 

Here, cr 2 is the projected pairwise peculiar velocity dispersion; the three-dimensional 
mean-square pairwise peculiar velocity (v 2 ) = 3<r 2 , since the mean pairwise peculiar velocity 
is taken to be zero. 

Given a value for the cosmological density parameter Qq, a may be estimated using the 
Cosmic Virial Theorem (Peebles 1980, Fisher et al. 1994b); 

_ 3H(z)Mz)QJr (zyr^ 
° { , ) - 46( 7 -l)(2- 7 )(4-7) ' (25) 
where Q relates the two- and three-point correlation functions, b is the linear bias factor 
and J depends only on 7 (J (7 = 1.7) = 4.14). For 7 close to 2, a is almost independent of 
separation, consistent with equation (|24]) . Equation (p5|) depends on the relation between 
the distributions of galaxies and matter through the (unknown) bias factor, and thus is of 
limited use as a probe of the true value of ft . Note that according to this model, a evolves 
as (1 + z)~ £ l 2 (holding 7 and Q constant). 
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3.3.2. Results 

The redshift space correlation function is calculated for the field sample using the 
redshift space analogue of equation 

RDD(s) . , 

where DD(s) and DR(s) are the number of data-data and data-random pairs with redshift 
space separations between s and s + As, respectively. The random catalog used here is the 
same as that used for computing the real space correlation function. Figure || shows £ (s) 
for the field sample, along with the predictions from equation (]23f) , using r = 2.1 h~ 1 Mpc, 
7 = 1.7 and a = 770 kins" 1 for Qq = 1, and r$ = 2.5 h~ 1 Mpc, 7 = 1.7 and a = 400 kms" 1 
for fl Q = 0.2 These values of a are computed using equation ( p5|) with Q = 6 = 1. Also 
shown are the curves given by equation (|23|) using a Gaussian pairwise peculiar velocity 
distribution with a = 140 kms" 1 ; this value corresponds to the mean velocity uncertainty 
in the sample of 100 kms -1 . 

As can be seen, the a = 770 kms -1 , Oo = 1 model overestimates the redshift space 
perturbations. A least-squares fit of equation (|23|) to the data, using equations (|24D and 
(|5|), with Q = 1, b = 1 and 7 = 1.7 yields r = 3.9tg;g /i -1 Mpc, inconsistent with the 
value derived from the projected correlation function data, tq = 2.l1i^ h~ 1 Mpc, with 
90% confidence. The a = 400 kms" 1 , Qo — 0.2 model matches the data more closely, 
consistent with the low favored by pavis fc Peebles (1983)| and [Fisher et al. (1994b)| . 



However, the observed is consistent with a model with zero pairwise peculiar velocity 
dispersion, and the redshift space distortions due solely to velocity measurement errors; 
more data are therefore needed for a precise determination of a. We conclude that the data 
are best modeled by a low density parameter; taking Q = 0.2, r = 2.5^0^ h~ 1 Mpc and 
o = 400 kms" 1 yields models which are consistent with both the observed w p (r p ) and 



4. Conclusions 

We have found that the physical correlation length for 0.21 < z < 0.53 is 
r = 2.li°J fr-^Mpc if tt = 1, implying e = 1.8±?;|. If fi = 0.2, r = 2.5^ h^Mpc and 
e = 0.8^{;3; the uncertainties are estimated using the bias-corrected bootstrap resampling 
method, with 300 resamplings. These results are consistent with earlier results obtained 
from angular surveys, which indicate rapid evolution (Efstathiou et al. 1991, Infante & 
Pritchet 1995). It is also consistent with the results from the CFRS (Le Fevre et al. 1996). 
This decrease in ro from its present value may be interpreted either as a real change in the 
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clustering of the observed galaxies, or as due to a weakly clustered population at z~0.36 
which is intrinsically faint at the present epoch. 

The projected pairwise peculiar velocity dispersion at z = 0.36, a = 770 km s -1 , 
predicted by the Cosmic Virial Theorem using Qq = 1 is inconsistent at the 90% confidence 
level with the observed redshift space correlation function. The Qq = 1 is therefore weakly 
rejected. The f2 = 0.2 prediction, a = 400 kms -1 , however, matches the data more closely. 
Thus, the relatively small redshift space distortions present favour low f2 as determined 
from the Cosmic Virial Theorem, consistent with the results of Pavis fc Peebles ( 1983 )| and 
Fisher et al. (1994b)| . 



More data is required to obtain a precise value for a, and larger scales need to be 
sampled in order to obtain a smooth redshift distribution, thereby removing uncertainties 
in £ due to density inhomogeneities on the scale of the sample diameter. A larger sample 
would also enable computations of the correlation function for subsamples based on galaxy 
colour or intrinsic brightness, which would help distinguish between the two possible sources 
of observed evolution described in section |3.2| . The CNOC2 redshift survey, presently 
in progress, will yield ~ 5000 high-accuracy redshifts in the range 0.15 < z < 0.7. This 
survey will contain enough objects, and sample sufficiently large scales, to permit accurate 
computations of the redshift space and real space correlation functions and their evolution 
at intermediate redshifts. 

We thank all participants of the CNOC cluster survey for assistance in obtaining and 
reducing these data. The Canadian Time Assignment Committee for the CFHT generously 
allocated substantial grants of observing time, and the CFHT organization provided the 
technical support which made these observations feasible. We gratefully acknowledge 
financial support from NSERC and NRC of Canada. 
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Fig. 1. — The MS1512+36 redshift sample. The angular scale has been expanded by a 
factor of~80. The objects within the central box are not included in the field sample. 

Fig. 2. — The redshift distribution ^-(z) for the redshift sample and corresponding random 
catalog. The solid line is ^ for the redshift sample; the distribution of the random catalog 
is indicated by the dotted line. The two vertical lines indicate the redshift range containing 
the 38 objects present in the redshift sample but not in the field sample. The random 
distribution has been normalized so as to have the same integral as the data distribution 
over the redshift range shown. 

Fig. 3. — The angular correlation function w{9) for the photometric and redshift samples. 
The values of w for the photometric sample are indicated by dots; the data for the redshift 
sample are indicated by crosses. The error bars are the 68.3% uncertainties estimated from 
300 bootstrap resamplings of the data. The solid line is given by the least-squares fit to the 
photometric data, with S in equation (9) fixed to 0.7. 

Fig. 4. — The counts in excess of random within a cylinder of radius 1 ft, _1 Mpc and length 
2ir, 2nn J p (l ft, _1 Mpc, 7r), for the field sample, for Q — 1 and 0.2. The error bars are the 
68.3% bootstrap confidence intervals. 

Fig. 5. — The projected correlation function w p {r p ) for the field sample, for Q = 1 and 0.2. 
The error bars are the 68.3% bootstrap confidence intervals. The solid line in each panel is 
from the least-squares fit of equation (17) to the data with S fixed to 0.7. 

Fig. 6.— The redshift space correlation function £(s) for the field sample, for Q = 1 
and 0.2. The error bars are the 68.3% bootstrap confidence intervals. The dashed line in 
each panel is the a power-law with 7 = 1.7 and r taken from the fit to the corresponding 
projected correlation function data. The solid curve in each panel is given by equation (23), 
with a = 770 km s -1 (Q = 1) and 400 km s -1 (Q = 0.2). The dotted curve is obtained 
using a Gaussian with a = 140 kms" 1 in place of the exponential in the velocity distribution 
model (24). 
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